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Abstract 

The availability of the double-differential charged-current neutrino cross section, measured by the MiniBooNE 
collaboration using a carbon target, allows for a systematic comparison of nuclear effects in quasi-elastic electron and 
neutrino scattering. The results of theoretical studies based on the impulse approximation scheme and state-of-the-art 
models of the nuclear spectral functions suggest that the electron cross section and the flux averaged neutrino cross 
sections corresponding to the same target and seemingly comparable kinematical conditions can not be described 
within the same theoretical scheme using the value of the nucleon axial mass obtained from deuterium measurements. 
I analyze the assumptions underlying the treatment of electron scattering data, and argue that the development of a 
new paradigm, suitable for application to processes in which the lepton kinematics is not fully determined, will be 
required. 
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1. Introduction 

The data set of Charged Current Quasi Elastic 
(CCQE) events recently released by the MiniBooNE 
collaboration [1| provides an unprecedented opportu- 
nity to carry out a systematic study of the double dif- 
ferential cross section of the process, 

V p + 12 C -> pT + X , (1) 

averaged over the neutrino flux. 

The charged current elastic neutrino-nucleon process 
is described in terms of three form factors. The proton 
(p) and neutron («) vector form factors, F^'"(Q 2 ) and 
F^'XQ 1 ) (Q 2 = -q 2 , q being the four-momentum trans- 
fer), have been precisely measured up to large values 
of Q 2 in electron-proton and electron-deuteron scatter- 
ing experiments, respectively (for a recent review, see, 
e.g., Ref.E). The (2 2 -dependence of the nucleon ax- 
ial form factor F A {Q 2 ), whose value at Q 2 — can be 
extracted from neutron /3-decay measurements, is gen- 
erally assumed to be of dipole form and parametrized in 



terms of the so called axial mass M A : 

F A (Q 2 )=g A (l + Q 2 /M 2 A y 2 . (2) 

The world average of the measured values of the ax- 
ial mass, mainly obtained from low statistics experi- 
ments carried out using deuterium targets, turns out to 
beM A = 1.03 ± 0.02 GeV EES, while the analyses 
performed by the K2K and MiniBooNE Q collab- 
orations using oxygen and carbon targets, respectively, 
yield M A ~ 1.2-^1 .35 GeV. 

It would be tempting to interpret the large value of 
M A reported by MiniBoonNE and K2K as an effective 
axial mass, modified by nuclear effects not included in 
the oversimplified Fermi gas model employed in data 
analysis. However, most existing models of nuclear 
effects (for recent reviews see Ref.[8]) fail to support 
this explanation, suggested by the authors of Ref.Q, a 
prominent exception being the model of Ref. 1 9 1 . 

Obviously, a fully quantitative description of the 
electron-scattering cross section, driven by the known 
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vector form factors, is a prerequisite for the under- 
standing of the axial vector contribution to the CCQE 
neutrino-nucleus cross section. 

Over the past two decades, the availability of a large 
body of experimental data has triggered the develop- 
ment of advanced theoretical descriptions of the nuclear 
electromagnetic response. The underlying scheme, 
based on nuclear many -body theory and realistic nuclear 
hamiltonians, relies on the premises that i) the lepton 
kinematics is fully determined and ii) the elementary in- 
teraction vertex can be extracted from measured proton 
and deuteron cross sections. 

The above paradigm has been successfully applied to 
explain the electron-nucleus cross section in a variety 
of kinematical regimes (for a recent review of the quasi- 
elastic sector see Ref. lflOll ). However, in view of the un- 
certainties associated with the energy of the incoming 
beam, the identification of the reaction mechanisms and 
the determination of the interaction vertex, its extension 
to the case of neutrino scattering may not be straightfor- 
ward. 

The comparison between theoretical calculations and 
data, presented in Section[2j suggests that the measured 
electron- and neutrino-nucleus cross sections can not be 
explained within the same scheme. The difficulties as- 
sociated with the description of the flux unfolded to- 
tal cross section are discussed in Section [3] Finally, 
in Section |4] I argue that the the paradigm succesfully 
employed to analyze electron-nucleus scattering is not 
suitable for application to processes in which the lepton 
kinematics is not fully determined, and must be signifi- 
cantly modified. 

2. Electron- and neutrino-nucleus interactions 

Electron-nucleus scattering cross sections are usually 
analyzed at fixed beam energy, E e , and electron scatter- 
ing angle, 6 e , as a function of the electron energy loss 
to. As an example, Fig. [TJshows the double differential 
cross section of the process 

e+ l2 C->e'+X, (3) 

at E e = 730 MeV and 6 e = 37°, measured at 
MIT-Bates |11|. The peak corresponding to quasi- 
elastic (QE) scattering, the bump at larger to, associated 
with excitation of the A-resonance, and the region in- 
between, where the cross section is mainly arising from 
processes involving meson exchange currents (MEC), 
are clearly recognizable. The three-momentum trans- 
fer |q| turns out to be nearly constant, its variation over 
the range shown in the figure being 5 5%. As a conse- 
quence, the cross section of Fig{T]can be readily related 



to the linear response of the target nucleus to a probe 
delivering momentum q and energy to, defined as 

S(q,to) =Y j \(n\Y j al +(i a k \0)\ 2 6(aj+E -E n ).(4) 

n k 

In the above equation, |0) and \n) represent the target 
initial and final states, with energies Eq and E„, respec- 
tively, while a k+(j and a^ are the nucleon creation and 
annihilation operators. 

The magnitude of the momentum transfer, 
|q| ~ 450 MeV, is large enough to make the im- 
pulse approximation (IA) scheme, in which the nuclear 
response of Eq.Q reduces to fT2l 

S M (q, to) = j d 3 k dEP,,(k, E)P p (k+q, oj-E) , (5) 

safely applicable fl3l . In Eq.(j5j, the hole and particle 
spectral functions, k , E ) and P p (k + q , cu — E), 
describe the energy and momentum distributions of the 
struck nucleon in the initial (hole) and final (particle) 
states, respectively. 

The solid line of Fig. [TJrepresents the results of a the- 
oretical calculation of the QE contribution lfl4l . carried 
out within the approach described in Refs. lfTBl [T6ll us- 
ing the hole spectral function of Ref.| 17 1 and the recent 
parametrization of the vector form factors of Ref.||4). 
Final state interactions (FSI) between the struck nucleon 
and the recoiling spectator system 1 15 1, whose main ef- 
fect is a ~ 10 MeV shift of the QE peak, have been also 
taken into account. 
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Figure 1: Inclusive electron-carbon cross section at beam energy 
E e = 730 MeV and electron scattering angle 8 e = 37°, plotted as a 
function of the energy loss ui 1141 . The data points are taken from 
Ref.fTTl. 

It is apparent that height, position and width of the 
QE peak, mostly driven by the energy and momentum 
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dependence of the hole spectral function, are well re- 
produced. 

Applying the same scheme employed to obtain the 
solid line of Fig. [JJto neutrino scattering leads to the re- 
sults shown in Fig. |2][ 14]. The data points represent the 
double differential CCQE cross section averaged over 
the MiniBooNE neutrino flux, whose mean energy is 
( E v ) - 788 MeV, plotted as a function of the kinetic 
energy of the outgoing muon, T^, at different values of 
the muon scattering angle 6^. The solid lines show the 
results (integrated over the cos 6^ bins) obtained using 
the same spectral functions and vector form factors em- 
ployed in the calculation of the electron scattering cross 
section of Fig. [JJ and a dipole parametrization of the 
axial form factor with M A = 1 .03 MeV. 

Comparison of Figs. [T] and [2] indicates that the 
electron and neutrino cross sections corresponding to 
the same target and seemingly comparable kinemati- 
cal conditions (the position of the QE peak in Fig. [TJ 
corresponds to kinetic energy of the scattered electron 
~ 610 MeV) can not be explained using the same the- 
oretical approach and the value of the axial mass re- 
sulting from deuterium measurements. In this instance, 
the paradigm of electron scattering appears to conspic- 
uously fail. 

The above conclusion, while being based on a calcu- 
lation carried out within the scheme of Refs. lT31[T6l . is 
largely model independent. Theoretical approaches pro- 
viding a quantitative description of the electron-nucleus 
cross section in the QE channel, are bound to pre- 
dict CCQE neutrino-nucleus cross sections significantly 
below the MiniBooNE data if the value of the axial 
mass is set to 1.03 GeV As a matter of fact, within 
the approach of Refs. lfT5l [TBI , the axial mass yelding 
the best ^ 2 -fit to the flux integrated MiniBooNE Q 2 - 
distribution, M A = 1.6 GeV, turns out to be even larger 
than that reported in Ref.Jl |. The authors of Ref.lfT4l 
have shown that this value of M A also explains the muon 
energy spectrum and angular distribution obtained from 
integration of the double differential cross seciton of 
Ref.|TJ. 

In Fig. [3] the theoretical results obtained using 
M A =1.6 GeV are compared to the distribution of the 
flux-averaged MiniBooNE events, plotted as a function 
of the reconstructed Q 2 , defined as 

Q 2 = 2E y E M cos 6^ - m\ , (6) 

where E M = T u + m^, = (E 2 - m 2 ) 1 ^ 2 and m^ is 
the muon mass. Note that the definition of Q 2 , besides 
the measured kinematical variables 7V and 0^, involves 
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Figure 2: Flux averaged double differential CCQE cross section mea- 
sured by the MiniBooNE collaboration [T], shown as a function of 
kinetic energy of the outgoing muon. The upper and lower panels 
correspond to to different values of the muon scattering angle. The- 
oretical results have been obtained using the same spectral functions 
and vector form factors employed in the calculation of the electron 
scattering cross section of Fig. ^ and a dipole parametrizaition of the 
axial form factor with M A = 1.03 MeV fT31 . 
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Figure 3: Flux averaged Q 2 -distribution obtained from the approach 
of Refs.(B][T6l using M A =1.6 GeV, compared to the data of Ref.| 1 1. 



the incoming neutrino energy E v , whose determination 
necessarily implies some assumptions. 
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3. Flux unfolded total cross section 

Figure [4] shows a comparison between the the Mini- 
BooNE ^mx unfolded total cross section and the results 
of the calculations of of Ref. [ 14 1 . It is apparent that 
in this case using M A = 1 .6 GeV leads to overestimat- 
ing the data in the region of high energy (E v > 800 
MeV), where the choice M A = 1 .35 GeV, resulting from 
he analysis of Ref.Hl, provides a better fit. The differ- 
ent pattern emerging from Fig. [4] compared to Fig. [3] 
clearly points to the uncertainty associated with the in- 
terpretation of flux averaged and flux unfolded data. 
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Figure 4: Flux unfolded total CCQE cross section, as a func- 
tion of neutrino energy. The dot-dash, solid and dashed lines 
have been obtained setting the value of the axial mass to 
Ma = 1.03, 1.35 and 1.6 GeV, respectively. The data are taken 
from Ref.Q] . 

A different scenario is suggested by the results of 
Ref. 15), whose authors obtain a quantitative account 
of the MiniBooNE flux unfolded total cross section 
without increasing M A . Within the model of Ref. fl9), 
the main mechanism responsible for the enhancement 
that brings the theoretical cross section into agreement 
with the data is multi-nucleon knock out, leading to n 
particle-n hole («p-«h) nuclear final states (n = 2, 3, . . .). 
Within the approach of Refs. |fT31[T6l , the occurrence of 
2p-2h final states is described by the continuum part of 
the spectral function, arising from nucleon-nucleon cor- 
relations ifTTl . It gives rise to the tail extending to large 
w, clearly visible in Fig. [T] However, its contribution 
turns out to be quite small (less than 10% of the inte- 
grated spectrum). The analysis of the momentum distri- 
bution sum rule indicates that the contributions of np-nh 
final states with n > 3 are negligibly small EOl. 

According to the philosophy outlined in this paper, in 
order to firmly establish the role of multi-nucleon knock 



out in CCQE neutrino interactions the model of Ref. [9 1 
should be thoroughly tested against electron scattering 
data. 

4. Summary and outlook 

The theoretical and experimental results discussed in 
this paper suggest that the main difference involved in 
the analysis of neutrino-nucleus scattering, as compared 
to electron-nucleus scattering, lies in the flux average. 

Unlike the electron cross section shown in Fig. [T] the 
flux averaged CCQE neutrino cross section at fixed en- 
ergy and scattering angle of the outgoing lepton picks up 
contributions from different kinematical regions, where 
different reaction mechanisms dominate. 

Consider, as an example, the cross section at 
cos^ = 0.75, the central value of the bin consid- 
ered in the upper panel of Fig. [2] and T M = 545 MeV, 
corresponding to the maximum of the measured cross 
section. Figure [5] displays the behavior of the Bjorken 
scaling variable x in these kinematical conditions and 
the MiniBooNe flux as a function of the neutrino energy 
E v . 
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Figure 5: Bjorken scaling variable x in the kinematical conditions 
corresponding to maximum of the cross section displayed in the upper 
panel of Fig. [2] plotted as a function of the neutrino energy E v . For 
comparison, the MiniBooNE neutrino flux is also shown. 

It turns out that x = 1, corresponding to quasi- 
elastic kinematics, and x = 0.5, corresponding to 
the dip region, are associated with neutrino energies 
E v - 788 MeV (the mean energy of the MiniBooNE 
flux) and 975 MeV, respectively, and that 

<t(£ v = 975 MeV) 



<6(£ v = 788 MeV) 



0.83 



(7) 
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It follows that the flux averaged cross section picks up 
the contributions of quasi-elastic scattering and MEC 
with about the same probability. Hence, it can not be de- 
scribed according to the paradigm successfully applied 
to electron scattering, based on the tenet that the lepton 
kinematics is fully determined. 

A new paradigm, suitable for studies of neutrino in- 
teractions, should be based on a more flexible model 
of nuclear effects, yielding a realistic description of the 
broad kinematical range associated with the relevant 
neutrino energies. 

Nuclear many-body theory provides a consistent 
framework for the development of such a model. Be- 
sides single- and multi-nucleon knock out, it should 
include the contributions of processes involving MEC, 
which are long known to provide a significant enhance- 
ment of the electromagnetic nuclear response in the 
transverse channel ETl . It has to be emphasized that 
quasielastic scattering and processes involving MEC 
lead to the same final state, and can not be distinguished 
in the MiniBooNE analysis. The occurrence of inelastic 
processes, leading to excitation of nucleon resonances 
and pion production should also be taken into account. 

As a final remark, it has to be emphasized that a 
great deal of information could be obtained applying 
the new paradigm to the analysis of inclusive neutrino- 
nucleus cross sections, preferably, although not neces- 
sarily, through direct implementation of the resulting 
nuclear model in the Monte Carlo simulation codes. 
This kind of analysis, which has been succesfully car- 
ried out for electron-nucleus scattering, may help to rec- 
oncile the different values of the axial mass obtained 
by different experiments, and would be largely unaf- 
fected by the problems associated with the possible 
misidentification of CCQE events, recently discussed in 
Ref. 11221 . 
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